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1.  Introduction 


Determination  of  the  average  flux  level,  or  the  thermally  radiant  intensity,  and  pulse 
length  from  an  aluminum  powder/  liquid  oxygen  Thermal  Radiation  Simulator  (TRS)  has 
generated  debate  because  of  the  nature  of  TRS  output  and  its  erratic  nature.  A  scheme 
is  needed  to  quantify  flux  against  time  output  into  single  values  for  TRS  characterization 
purposes. 

Past  attempts  to  shape  the  thermal  pulse  from  a  TRS  to  match  that  of  a  nuclear  weapon 
have  not  been  entirely  successful.  The  TRS  rectangular  pulse  is  the  most  widely  used  and 
accepted.  Current  technology  permits  such  pulses  to  be  produced  with  rapid  rise  and  fall 
times  (less  than  100  msec,  to  full  flux  capacity). 

The  problem  is  how  to  quantify  the  thermal  output  into  simple  values  for  comparisons 
against  other  TRS  performances  or  for  predicting  target  response.  The  TRS  can  be  pro¬ 
grammed  to  turn  on  and  off,  but  the  thermal  output  is  still  a  widely  variant  distortion  of  a 
rectangular  pulse.  Depending  on  the  location  of  the  heat  sensing  device,  the  variation  can 
exceed  up  to  20%  of  the  observed  mean  flux.  Determining  values  for  average  flux  and  pulse 
time  are  difficult  to  establish  with  specificity. 

Four  schemes,  for  determining  the  average  flux  and  pulse  duration  are  presented.  They 
are  the  1)  Full  Width  Method,  2)  Full  Width  at  Half  Maximum  Method,  3)  Moments  Match¬ 
ing  Method  and  4)  the  Fourier  Averaging  Method.  Each  method  is  presented  with  its  respec¬ 
tive  advantages  and  limitations.  The  results  are  then  compared  to  the  predicted  analytical 
and  numerical  responses  of  different  targets.  These  methods  all  represent  attempts  to  quan¬ 
tify  the  TRS  data  in  a  consistent  manner  for  comparison  with  other  TRS  data. 


2.  Thermal  Radiation  Simulators 


There  exist  at  least  four  different  methods  currently  used  for  nuclear  thermal  radiation 
simulation  [1].  Solar  furnaces,  such  as  the  White  Sands  Missile  Range,  New  Mexico,  and 
Centre  d’Essais  d’Odeillo,  France,  consist  of  a  large  mirror  that  tracks  the  relative  movement 
of  the  sun  with  the  earth  and  reflects  the  suns  rays  onto  a  parabolic  array  of  concentrator 
mirrors  .  The  sunlight  is  focused  into  a  small  area  of  irradiation.  Because  of  the  high 
flux  levels  achievable,  and  the  nature  of  the  system,  nuclear  event  pulses  may  be  imitated 
with  accuracy.  Unfortunately,  the  concentrator  only  irradiates  a  small  area,  and  the  facility 
cannot  be  used  for  large  targets. 

The  Xenon  Flashlamp  Facility  and  Quartz  Lamp  Bank  Facilities  at  Tri-Service  Thermal 
Radiation  Test  Facility,  Wright-Patterson  Air  Force  Base,  Ohio,  have  consistent  outputs,  but 
are  limited  by  the  amount  of  target  area  that  can  be  exposed.  The  Xenon  Flashlamp  can 
simulate  the  nuclear  pulse,  but  like  the  solar  facility,  can  only  expose  a  target  4  by  4  in.  [lj. 
The  Quartz  Lamp  Bank  is  similar  to  the  TRS  in  that  its  output  is  also  a  square  wave.  It  is 
plagued  by  a  small  irradiation  surface,  a  relatively  low  flux  output  level  (10  inch  long  by  6 
inch  high,  55  cal/cm2s)  and  large  electrical  energy  requirements. 
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The  aluminum/liquid  oxygen  Thermal  Radiation  Simulator  uses  the  latent  energy  re¬ 
leased  by  radiation  of  the  condensing  particles  of  aluminum  oxide  [3].  The  aluminum  powder 
is  jetted  into  the  air  at  high  velocity  along  with  saturated  liquid  oxygen.  The  mixture  passes 
an  ignition  point  and  combusts.  This  process  creates  a  large  wall  of  flame  radiating  at  a  high 
temperature.  Unfortunately,  the  TRS  flame  is  extremely  sensitive  to  boundary  conditions, 
yielding  variation  in  the  output.  As  shown  in  Figure  1,  the  TRS  output  does  yield  some 
average  flux  over  some  definite  time  frame. 


Figure  1.  Typical  TRS  output  of  the  BRL  Facility,  APG,  MD 


3.  Data  Analysis  Schemes 


The  output  of  the  Thermal  Radiation  Simulator  is  difficult  to  control.  Regulation  of  the 
aluminum  powder  flow  to  the  combustion  nozzle  is  difficult  without  a  sophisticated  control 
feedback  system.  The  BRL  TRS  procedure  is  to  preset  the  aluminum  containment  vessel 
pressures.  A  valve  at  the  bottom  of  the  vessel  opens  and  closes  aluminum  flow  to  the  system. 
A  valve  at  the  combustion  nozzle  controls  the  direction  of  aluminum  flow.  Aluminum  flows 
either  to  a  waste  recovery  system  or  to  the  combustion  chamber.  During  TRS  operation, 
the  aluminum  powder  is  allowed  to  flow  in  the  waste  recovery  system  for  a  few  seconds,  and 
then  diverted  into  the  combustion  chamber  with  minimal  disturbance  to  the  steady  flow. 
The  result  is  a  more  controlled  rectangular  thermal  pulse.  While  the  flux  record  does  not 
match  that  of  a  nuclear  pulse,  the  fluence  is  matched,  effectively  exposing  a  target  to  an 
equivalent  thermal  environment. 

Since  the  TRS  output  is  rectangular,  several  schemes  have  been  developed  to  quantify 
the  output  as  an  ideal  rectangular  thermal  pulse  with  defined  amplitude  and  duration.  The 
amplitude  can  be  described  as  the  average  flux.  The  duration  is  the  time  of  thermal  output 
at  the  defined  average  flux,  or  pulse  width. 


3.1  Full  Width  Method 

The  Full  Width  method  defines  the  pulse  width  of  the  flux  record  as  the  time  from  the 
initiation  of  the  TRS,  until  the  flux  level  has  dropped  below  some  minimum  point  close  to 
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a)  ' 

Figure  2.  Comparisons  between  a)Flux  record  and  b)  ideal  Rectangular  pulse 

the  baseline.  The  amplitude  is  then  determined  by  dividing  the  fluence,  or  the  integrated 
flux  record,  by  the  pulse  width.  The  beginning  of  a  flux  record  ordinarily  is  well  defined, 
but  the  termination  point  is  vague  due  to  hot  cloud  residue  still  present  to  the  target  well 
after  shut-down.  This  leaves  a  decaying  trail  on  the  calorimeter  data  and  prolongs  the 
termination  point.  As  seen  in  Figure  3,  the  determination  of  the  pulse  width  usually  leads 
to  low  amplitude  levels,  poorly  characterizing  the  TRS  data. 


Figure  3.  Comparisons  between  TRS  Data  and  Full  Width  Method 


3.2  Full  Width  at  Half  Maximum  Method 

This  method  is  very  much  a  corrected  Full  Width  method.  The  Full  Width  method 
is  modified  by  redefining  the  pulse  width  onset  and  termination  points.  These  points  are 
defined  as  where  the  flux  record  initially  reaches  more  than  one  half  the  maximum  value 
and  where  the  flux  level  drops  below  this  value  again.  This  method  improves  the  ideal 
characterization  because  it  is  not  prolonged  by  the  residual  burn  off. 

Large  variations  in  TRS  output  affect  the  repeatability  of  the  TRS  maximum.  Deter¬ 
mination  of  the  half  maximum  value  changes  between  records.  Because  of  the  amount  of 
this  variation,  crossover  to  below  half  maximum  could  occur  more  than  once.  This  presents 
a  problem  in  some  algorithms  for  establishing  the  parameters  but  can  be  overcome.  A  draw- 
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back  to  this  method  is  the  exclusion  of  data  outside  the  defined  pulse  length  resulting  in  a 
low  fluence  calculation. 


Figure  4.  Comparisons  between  TRS  Data  and  Full  Width  Method  at  Half  Maximum 


3.3  Moments  Matching  Method 


The  Moment  Matching  Method  can  establish  the  parameters  of  the  rectangular  pulse 
from  an  arbitrary  pulse  [2]  in  terms  of  its  zeroth,  first  and  second  moments.  The  n’th  moment 
is  defined  as 

/+oo 

f{x)tndt  (1) 

-OO 


fiThe  rectangle  parameterized  is  the  one  which  most  closely  matches  the  same  moments. 
The  amplitude  A  and  the  pulse  width  pw  can  be  found  centered  about  a  point  in  time,  tc, 
from  the  following 


Mi 


tc~  M, 


pw  =  2\l3{Mi-t'2) 


A  = ^ 
pw 


(2) 

(3) 

(4) 


The  zeroth  moment, M0,  is  the  fluence  of  the  record,  and  establishing  pw  in  somewhat  rigorous 
fashion.  The  limitation  of  this  characterization  scheme  is  that  pw  is  dependent  on  the  second 
moment,  Af2.  The  second  moment  is  extremely  sensitive  to  the  range  of  data  being  observed. 
As  the  range  of  data  points  is  increased,  the  instrumentation  fine  noise  will  quickly  affect 
the  second  moment  value  because  of  the  t 2  term.  The  consequence  is  uncertainty  regarding 
the  calculated  pulse  width.  The  TRS  data  set  must  be  carefully  bounded  to  eliminate  fine 
noise  outside  what  is  considered  pertinent  data,  before  analysis  can  be  performed. 


3.4  Fourier  Averaging  Method 

The  Fourier  Averaging  Method  is  a  completely  formulated  method  of  characterizing 
a  TRS  data  set,  requiring  no  estimation  or  handwork.  The  idea  is  much  the  same  as  the 
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Moments  Matching  method  in  that  the  data  set  is  subjected  to  a  process,  then  compared  to 
a  rectangular  pulse  subjected  to  the  same  process.  The  rectangular  pulse  that  most  closely 
matches  that  data  defines  the  characterization  parameters. 

The  Fourier  transform  of  a  time  dependent  equation  into  the  frequency  domain  is  defined 
as  such 

/+°°  , 

F{t)-"’'it  (6) 

■oc 

When  the  integration  is  performed  at  /  =  0,  the  result  is  the  fluence.  The  method  for 
finding  the  pulse  width  is  done  by  finding  the  first  point  in  the  frequency  domain  where 
the  magnitude  of  the  transform  is  a  minimum.  This  frequency  point  is  equivalent  to  the 
reciprocal  of  the  pulse  width.  For  a  true  rectangular  pulse  [4] 

PW  =  JQ  (7) 

Fuo)  =  0  (8) 

The  Fourier  transform  of  a  TRS  data  set  will  yield  approximately  the  same  transform. 


Fourier  Transform 


Figure  5.  Fourier  Transform  of  a  Rectangular  Pulse 

By  simply  finding  the  first  minimum  point  of  the  transformed  record,  the  pulse  width  is 
established.  The  TRS  data  set  is  then  fully  characterized.  The  transform  is  relatively 
insensitive  to  the  range  of  data  used.  Because  of  the  insensitivity,  data  points  can  be  bypassed 
to  increase  the  speed  of  characterization. 

The  Fourier  transform  will  break  into  two  parts,  real  and  imaginary.  The  magnitude  is 
the  root  of  the  sum  of  the  squares  of  the  two  parts.  At  zero  hertz,  the  energy  density  will 
be  at  a  maximum.  This  is  also  the  fluence  of  the  TRS  record.  As  the  frequency  increases, 
the  resulting  integration  will  decrease  until  eventually  a  minimum  is  reached.  Since  the 
magnitude  is  the  root  of  the  sum  of  the  squares,  its  value  will  never  be  negative.  For  a 
true  rectangular  pulse,  this  first  minimum  will  be  zero.  Since  TRS  is  a  distortion  of  the 
rectangle,  the  minimum  may  not  be  zero,  but  the  minimum  is  still  clearly  defined.  The  pulse 
width  is  defined  as  the  reciprocal  of  the  frequency  where  the  magnitude  of  the  transform  is 
a  minimum.  For  a  more  detailed  treatment,  see  appendix  A. 

Finding  the  minimum  is  done  by  iteration,  as  shown  in  appendix  C.  Since  the  minimum 
is  so  clearly  defined,  few  iterations  are  required.  The  process  can  be  made  more  rapid  by  first 
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estimating  the  pulse  width  and  concentrating  the  method  in  that  area,  but  is  not  necessary. 
The  result  is  a  method  of  characterizing  the  TRS  thermal  pulse  without  any  estimations  or 
handwork.  The  technique  is  completely  rigid,  and  the  result  repeatable  with  precision. 


4.  Simulation  of  a  Target  Response 


To  evaluate  how  accurate  each  method  is  in  describing  the  TRS  output,  a  computer 
generated  thermal  response  simulation  was  performed.  An  imaginary  aluminum  plate  was 
exposed  to  actual  TRS  data  generated  at  BRL,  and  a  temperature  profile  calculated.  The 
four  characterization  methods  were  used  to  derive  idealized  thermal  pulses.  The  aluminum 
plate  was  exposed  to  the  thermal  pulses.  The  resulting  temperature  profiles  were  compared 
with  the  original  profiles. 

The  numerical  method  was  shown  to  work,  and  to  what  degree  of  accuracy  by  comparing 
the  numerical  results  to  an  analytical  solution  of  the  same  problem.  Once  satisfied  the 
numerical  method  worked,  the  comparisons  of  the  temperature  profiles  where  performed. 


4.1  Analytical  Solution 

To  evaluate  the  validity  of  the  characterization  schemes,  a  numerical  simulation  of  the 
TRS  data  and  the  results  of  the  four  schemes  are  compared.  To  demonstrate  the  validity  of 
the  numerical  method,  an  analytical  solution  of  the  same  problem  is  solved.  The  results  of 
the  analytical  solution  can  be  compared  to  those  of  the  numerical  method.  The  simulation 
target  is  an  infinite  slab  of  aluminum  with  defined  thickness.  To  simplify  the  calculations, 
the  following  assumptions  are  made; 

All  thermal  radiant  energy  is  absorbed. 

Re-radiation  from  the  slab  is  neglected. 

All  material  properties  are  assumed  constant  regardless  of  temperature. 

No  other  heat  transfer  processes,  such  as  convection,  are  considered. 

The  slab  is  infinitely  large,  therefore  limiting  the  calculation  to  one  dimension  (no  edge 
effects). 


The  following  is  the  transient  heat  conduction  equation  with  the  listed  symbols  indicated  in 
Figure  6. 

L  =  slab  thickness 
0  =  temperature 

x  =  distance  from  adiabatic  surface 
t  =  time 
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Figure  6.  Infinite  slab  considered  for  target  response 

k  =  thermal  conductance 
Cp  =  specific  heat 
p  =  material  density 
q  =  radiant  flux  per  unit  area 
a  =  thermal  diffusivity, 

The  heat  conduction  equation  is 

d2Q  IdQ 
dx 2  a  dt 

With  the  following  boundary  and  initial  conditions, 


1) t  =  O,0(X>O  =  O 

2) x  =  0,f  =0 

3) x  =  i,f  =  f 

the  solution  to  the  heat  conduction  equation  is  given  by  [5] 

+  uo) 

This  represents  the  first  part  of  the  solution  for  target  response,  exposing  the  target  to  the 
thermal  pulse.  The  second  part  of  the  solution  is  the  target  response  after  thermal  exposure. 
This  is  done  by  solving  the  heat  conduction  equation  for  the  slab  with  some  arbitrary  initial 
temperature  distribution  throughout  the  slab,  f(x).  Given  the  following  initial  conditions 
and  boundary  conditions,  and  selecting  a  new  time  variable  to  take  into  account  the  offset 
of  the  post-pulse  time, 

w  =  t  —  pw  (11) 

With  the  following  boundary  and  initial  conditions, 

l)tn  —  0,  Q{x,w)  ~  f(x) 

f(x)  is  the  initial  temperature  distribution 


7 


in  the  slab,  i.  e.  Equation  10  evaluated  at  t  =  pw 


2)*  =  0,H  =0 

3  )*  =  £,§?  =  -* 

the  solution  to  the  heat  conduction  equation,  again  by  [5],  is 

ge-on  Sir’t/L’C05( 

n=l 

For  a  more  detailed  explanation  about  how  this  solution  was  determined,  refer  to  appendix  C. 
This  analysis  is  good  for  the  transient  temperature  distribution  after  the  target  is  exposed  to 
the  TRS  pulse.  One  more  temperature  that  can  be  evaluated  is  the  final  uniform  temperature 
distribution  in  the  target  well  after  the  thermal  pulse.  The  result  can  be  easily  derived  as 
follows: 

A  -SJ^L  (]*) 

equilibrium  r  K*-0/ 

pCpL 

With  the  analytical  solutions  available,  the  numerical  method  can  be  validated. 


nnx  [L  no. 

—  )J  f(M)Cos(—)ds  (12) 


4.2  Numerical  Solution 


The  numerical  solution  is  used  to  observe  the  effect  of  actual  TRS  data  on  the  target. 
An  explicit  formulation  of  the  difference  equation  is  used  to  solve  the  transient  tempera¬ 
ture  history  given  the  same  assumptions  and  boundary  conditions  as  used  in  obtaining  the 
analytical  solution.  The  formulation  is  as  follows; 


T+  =  T,  +  q 


—  \Si 

Ax  L  k 


+  £ 


Tj-T, 

Ax 


(14) 


Ti  is  the  temperature  of  node  i  in  the  plate.  T+  is  the  explicit  prediction  of  the  temperature 
after  some  time  increment  At ,  and  Ai  is  the  incremental  node  thickness. 

The  TRS  data  was  stored  in  an  ASCII  file.  A  BASIC  program  was  written  and  run  on 
a  personal  computer  to  perform  the  explicit  calculation. 

The  simulation  used  in  this  report  was  a  15  centimeter  thick  aluminum  slab.  The  same 
assumptions  were  applied  as  in  the  analytical  solution.  It  was  felt  that  this  type  of  simulation 
would  be  most  appropriate  since  target  response  of  thermally  thick  targets  is  more  affected 
by  flux  levels  than  thermally  thin  targets  [7].  The  condition  for  thermal  thickness  occurs 
when 

a  ^<0.2  (15) 

where  a  is  the  thermal  diffusivity,  tmax  is  the  time  at  which  maximum  irradiation  from 
a  nuclear  event  occurs,  and  L  is  the  plate  thickness.  Thermally  thin  plates  are  affected 
mainly  by  the  fluence  of  the  source  irradiator.  If  a  plate  is  thermally  thin,  the  temperature 
distribution  is  essentially  uniform,  and  the  problem  is  treated  as  a  lumped  capacitance. 
Targets  exposed  to  different  events  with  equal  fluence  will  react  equally.  Thermally  thick 
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Table  1.  Results  of  numerical  and  analytical  comparisons 


Peak  Temp.  Front  Temp  of  Slab 


A 

pw 

Num. 

Anal. 

6 

Num. 

Anal. 

6 

40 

1 

386.9 

384.9 

2.0 

305.9 

304.6 

1.3 

50 

0.5 

374.2 

375.1 

-0.9 

303.5 

302.0 

-0.6 

60 

5 

604.8 

584.9 

19.9 

342.4 

334.5 

7.9 

30 

10 

517.8 

501.5 

16.3 

342.3 

334.5 

7.8 

40 

5 

503.2 

489.9 

13.3 

328.3 

323.0 

6.3 

100 

3 

689.3 

667.8 

21.5 

342.4 

334.5 

7.9 

100 

0.5 

448.5 

450.2 

-1.7 

307.1 

305.8 

1.3 

50 

1 

408.6 

406.2 

2.4 

307.2 

305.8 

1.4 

25 

2 

378.6 

375.1 

3.5 

307.1 

305.8 

1.3 

20 

5 

401.6 

395.0 

6.6 

314.2 

311.5 

2.7 

targets  respond  differently  to  events  with  different  flux  levels,  even  though  the  fluence  may 
be  equal.  These  targets  experience  higher  temperatures  on  the  exposed  surface  for  higher 
flux  events.  For  equal  fluence  exposures,  long  after  the  thermal  event,  the  target  temperature 
will  distribute  to  an  equilibrium  level. 

The  results  of  several  different  rectangular  thermal  pulses  were  used  to  check  the  numer¬ 
ical  results  against  the  analytical.  The  two  methods  yielded  the  same  results  consistently, 
indicating  the  numerical  method  worked  reasonably  well.  A  computer  simulation  of  the 
described  target  was  exposed  to  ten  random  thermal  pulses.  The  resulting  temperature  pro¬ 
files  yielded  a  peak  temperature  and  a  uniform  temperature  after  a  long  time.  These  two 
temperatures  are  compared  in  Table  1.  The  small  6  is  the  difference  between  the  numerical 
and  analytical  results. 

Comparisons  are  poor  for  long  duration  burns  due  to  error  growth  of  the  finite  difference. 
The  TRS  data  used  in  this  analysis  are  all  short  duration,  within  two  seconds.  From  Table 
1,  equivalent  comparisons  are  quite  close,  within  1%.  The  magnitude  of  the  flux  level  also 
amplifies  the  numerical  error,  but  the  data  used  in  the  evaluation  are  well  below  fifty  calories, 
again  maintaining  a  good  affin;  y  between  the  methods. 
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4.3  Target  and  Thermal  Radiation  Simulator  Data 


The  target  used  in  this  simulation  is  an  aluminum  slab  with  the  following  material 
properties. 

L  =  thickness,  15  centimeters 

p  -  density,  2,707  $ 

Cp  =  specific  heat,  0.896  -j—R 

k  =  thermal  conductance,  204  -",7; 

a  =  thermal  diffusivity,  8.418x10-5  “ 

The  assumptions  are  the  same  as  for  the  analytical  solution. 

All  thermal  radiant  energy  is  absorbed. 

Re-radiation  from  the  slab  is  neglected. 

All  material  properties  are  assumed  constant  regardless  of  temperature. 

No  other  heat  transfer  processes,  such  as  convection,  are  considered. 

The  slab  is  infinitely  large,  therefore  limiting  the  calculation  to  one  dimension  (no  edge 
effects). 

The  TRS  data  was  obtained  from  the  U.S.  Army  Ballistic  Research  Laboratorys  alu¬ 
minum/liquid  oxygen  Thermal  Radiation  Simulator.  Five  such  records  are  presented  ranging 
from  good  rectangular  pulses  to  poor  representations.  This  is  to  show  the  extremes  that  can 
be  found  in  real  life,  and  compare  the  results  of  the  four  characterization  schemes.  Figure  7 
shows  the  five  TRS  record  traces.  The  five  data  records  represent  the  thermal  output  sensed 
by  a  Medtherm  Calorimeter  placed  100  cm.  from  the  nozzle  center  and  29  cm.  above  the 
nozzle  top,  as  depicted  in  Figure  8.  This  is  the  reference  station  used  in  evaluating  and 
comparing  the  TRS  records  at  BRL.  The  incongruity  in  output  is  caused  by  certain  changes 
experimented  with  in  an  effort  to  improve  the  consistency  and  efficiency  of  TRS  burns. 

Trace  1  and  3  are  the  worst  case  data  sets.  Both  are  considered  poor  in  quality,  but 
are  presented  here  to  show  the  difficulties  in  achieving  consistent  parametric  levels.  Trace 
2  is  an  anomalous  TRS  burn  with  large  variants,  and  Traces  4  and  5  are  considered  close 
approximations  of  a  rectangular  pulse. 
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Figure  7.  Five  U.S.  BRL  TRS  Record  Traces  used  for  Characterization  scheme  comparisons. 


Figure  8.  Spatial  location  method  used  in  locating  calorimeter  position  relative  to  BRL  TRS 
Nozzle. 


11 


5.  Results  of  the  Parameterization  of  the  Five  TRS  Data 

Traces 


A  Fortran  program  was  written  to  evaluate  the  characteristic  parameters  concurrently. 
All  four  methods  were  evaluated  using  the  same  data  range  to  maintain  equal  fluence,  except 
in  the  Fuil  Width  Half  Maximum  method,  by  definition.  The  results  of  establishing  the  TRS 
data  parameters  is  listed  in  Table  2.  The  parameters  listed  below  were  used  to  generate 


Table  2.  TRS  Data  Parameters  as  established  by  the  Characterization  Schemes. 

Full  Width  Half  Max.  Mom.  Match.  Fourier  Ave. 


Flux 

Time 

Flux 

Time 

Flux 

Time 

Flux 

Time 

Trace  1 

28.01 

0.965 

31.02 

0.810 

31.39 

0.861 

32.03 

0.844 

Trace  2 

25.90 

1.275 

36.38 

0.835 

36.47 

0.905 

38.90 

0.849 

Trace  3 

15.23 

0.700 

19.75 

0.465 

18.50 

0.576 

19.87 

0.537 

Trace  4 

32.96 

1.800 

38.05 

0.89 

38.82 

0.917 

39.80 

0.895 

Trace  5 

34.61 

1.285 

38.85 

uoo 

39.93 

1.114 

40.67 

1.094 

twenty  ideal  simulation  rectangular  thermal  pulses,  and  by  the  same  method  as  before,  used 
to  derive  by  explicit  numerical  calculation,  the  temperature  profile  of  the  aluminum  target. 
These  profiles  were  then  compared  to  the  simulation  that  used  real  TRS  data.  Figures  9,  10 
and  11  shows  the  comparisons.  The  results  of  the  comparisons  show  the  difference  between 
the  ideal  pulse  and  the  TRS  pulse  is  small,  except  for  the  Full  Width  method.  The  consistent 
”best  fits”  are  the  Moments  Matching  and  Fourier  Average  methods.  The  Full  Width  method 
consistently  underestimates  the  flux  level,  resulting  in  a  lower  surface  temperature.  The  Full 
Width  at  Half  Maximum  consistently  under  calculates  the  fluence,  resulting  again  in  lower 
surface  temperature  rise  as  well  as  a  lower  final  temperature  distribution. 

The  advantage  of  the  Moments  Matching  is  the  calculation  for  the  parameters  is  rel¬ 
atively  straight  forward.  Only  three  integrations  are  necessary,  and  the  parameters  are 
computed  from  the  resulting  moments.  The  problem  is  that  by  changing  the  range  of  data, 
or  the  data  set  size,  used  in  the  computation,  the  parameters  change.  The  second  moment 
is  extremely  sensitive  to  line  noise,  so  care  must  be  taken  to  filter  the  data. 

The  Fourier  Average  method  is  insensitive  to  noise.  Figure  12  shows  the  effect  of  using 
the  entire  data  record  for  Trace  4  compared  to  the  carefully  ranged  data  set.  The  Full 
Width  method  result  in  a  pulse  length  of  the  entire  record  of  low  average  flux.  The  results 
is  a  different  temperature  profile  as  shown  below.  The  Moments  Matching  experiences  a 
pulse  lengthening  since,  by  Equation  3  pulse  length  is  computed  using  the  second  moment, 
as  well  as  a  slight  increase  in  the  fluence.  The  Fourie”  Average  shows  only  a  slight  increase 
in  the  flux  level  due  to  summing  the  noise  into  the  fluence,  but  the  pulse  length  essentially 
remains  the  same.  In  effect  the  Fourier  Average,  method  establishes  the  same  pulse  width, 
despite  the  increase  in  data  record.  The  drawback  to  using  the  Fourier  Average  method 
is  that  the  computation  involves  an  iterative  process  using  integrations.  This  will  result  in 
longer  computations  on  the  computer.  Depending  on  the  size  of  the  data  set,  that  time  may 
be  insignificant. 
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Since  the  Fourier  Average  method  is  insensitive  to  high  frequency  fluctuations,  one 
could  select  a  fraction  of  the  data  (pick  one  out  of  every  ten  sample  points)  to  reduce  the 
computation  time  for  pulse  length.  Figure  13  shows  the  comparison  of  the  Fourier  transform 
of  Trace  2  and  Trace  3.  The  original  sample  rate  was  200  points  per  second.  Data  decimation 
of  ten  to  one  reduces  the  data  sample  rate  to  20  points  per  second.  Trace  3  showed  a  wider 
variation  in  the  derivative  of  the  Fourier  transform  of  the  magnitude,  but  the  end  result 
was  an  extremely  close  match  for  average  flux  and  pulse  width.  One  must  add  that  the 
same  decimated  data  yield  equally  close  results  when  used  with  the  other  characterization 
schemes,  as  demonstrated  below  in  Table  3. 


Table  3.  Comparisons  of  Characterization  Schemes  results  between  full  data  set  and  deci¬ 
mated  data  set. 


Flux  (cal / cm2 sec)  Pulse  (sec) 


Full 

Decimated 

Full 

Decimated 

Full  Width 

Trace  2 

25.9 

25.43 

1.28 

1.3 

Trace  3 

15.23 

14.33 

Half  Maximum 

Trace  2 

36.38 

36.85 

0.835 

0.925 

Trace  3 

19.75 

19.65 

0.47 

0.5 

Moments  Matching 

Trace  2 

36.47 

35.98 

0.90 

0.919 

Trace  3 

18.5 

18.32 

0.576 

0.587 

Fourier  Average 

Trace  2 

38.9 

38.38 

0.849 

0.861 

Trace  3 

19.87 

19.73 

0.537 

0.545 
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Tempena+une  Profile 


-  Trace  1  —  Full  Uidth 

—  Half  Max i mum  -  Moments  Match i ng 

-  -  Fourier  fiverag 


Temperature  Profile 


-  Trace  2  —  Full  Uidth 

—  Half  Max i mum  —  Moments  Match i ng 

-  -  Fourier  Pverao 


Figure  9.  Comparisons  of  the  results  of  numerically  computed  temperature  profiles  on  the 
target  surface  using  Trace  1  and  Trace  2  with  their  respective  idealized  thermal  pulses. 
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Temperature  Prof i le 


- Trace  3  ---rail  Width 

—  Half  Max i mum  —  Moments  Match i ng 

-  -  Fourier  Av/erao 


Temperature  Profile 


-  Trace  4  -  Ful I  Width 

—  Half  Maximum  -  Moments  Matching 

■  -  Fourier  Averao 


Figure  10.  Comparisons  of  the  results  of  numerically  computed  temperature  profiles  on  the 
target  surface  using  Trace  3  and  Trace  4  with  their  respective  idealized  thermal  pulses. 
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Tempenatune  Prof i le 


-  Trace  5  - Ful  I  Width 

—  Half  Maximum  —  Moments  Matching 

-  -  Fourier  ftverao 


Figure  11.  Comparisons  of  the  results  of  numerically  computed  temperature  profiles  on  the 
target  surface  using  Trace  5  with  its  respective  idealized  thermal  pulses. 
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Temperature  Profile 


-  Trace  4 

—  Half  Maximum 
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Full  Uidth 
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a) 


Temperature  Profile 


-  Trace  4  —  Full  Width 

—  Half  Maximum  —  Moments  Matching 

-  -  Fourier  Averoo 


*) 

Figure  12.  Comparisons  of  the  results  of  numerically  computed  temperature  profiles  on  the 
target  surface  using  Trace  4  and  the  respective  four  idealized  thermal  pulse  taken  over  a)  a 
selected  range  of  data  and  b)  the  entire  range  of  data. 


17 


Fourier  Transform 


o  .5  i  i.5  a  a.s  3 


Frequency,  Hz 
Trace  2 

—  Magnitude  red.  -  d(aag>/df  red.  —  Magnitude 

- d<»ag)/df 


a) 


Fourier  Transform 


Frequency,  Hz 
Trace  3 


- Magnitude  red.  —  d<mog>/df  red.  —  Magnitude 

--  d(»og)/df 

b) 


Figure  13.  Comparisons  of  the  Fourier  Transform  of  a)  Trace  2  and  decimated  Trace  2,  b) 
Trace  3  and  decimated  Trace  3. 
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6.  Conclusions 


After  evaluating  the  four  parameterization  schemes  against  a  computer  thermal  response 
target  simulation,  the  Fourier  Average  technique  emerges  as  the  most  consistent  method. 
Regardless  the  size  of  the  data  set  and  noise,  the  parameters  of  the  TRS  data  are  calculated 
with  greater  consistency  than  with  any  other  method.  Decimating  the  data  had  little  effect 
on  the  output.  The  simulated  target  temperature  profiles  of  the  calculated  ideal  thermal  flux 
records  showed  the  Fourier  parameterized  pulse  was  consistent  with  the  actual  TRS  data 
simulation. 

With  the  increase  of  data,  the  pulse  width  of  the  Fourier  Average  calculated  ideal  pulse 
showed  a  slight  increase  in  average  flux  level.  The  Moments  Matching  methods  showed  a 
large  deviation  of  pulse  length  from  the  original,  carefully  isolated  data  set.  This  is  because 
the  second  moment  is  particularly  sensitive  to  instrumentation  line  noise.  An  inflated  value 
for  the  second  moment  will  result  in  an  increase  of  the  pulse  length.  The  Other  methods 
remained  the  same  by  definition.  The  Full  Width  and  Full  Width  at  Half  Maximum  methods 
ignores  all  data  outside  the  selected  data  range.  The  slight  increase  in  average  flux  was  due 
to  the  line  noise  is  also  being  integrated.  The  effect  is  reflected  in  the  Moments  Matching 
method  as  well. 

The  drawback  of  the  Fourier  Average  technique  is  the  iterative  method  used  to  find 
the  pulse  width.  With  an  extremely  large  data  set,  this  could  result  in  long  computation 
time.  Elimination  of  ninety  percent  of  the  data  set  did  not  affect  the  method  results  by 
more  than  1  %.  This  would  increase  the  speed  of  computation  with  little  to  no  effect  on  the 
outcome.  A  typical  Thermal  Radiation  Simulator  data  set  consists  of  about  200  points  per 
second.  For  a  three  second  data  set,  the  calculation  time  difference  on  a  personal  computer 
is  indiscernible. 

The  result  of  the  Fourier  Average  generated  ideal  thermal  pulse  against  the  actual 
TRS  data  in  the  aluminum  plate  simulation  showed  good  consistency.  The  Fourier  Average 
method  was  not  consistently  better  at  parameterizing  an  ideal  match.  For  different  TRS 
data  sets,  some  other  methods  worked  better.  This  is  due  to  uniqueness  in  a  data  set.  The 
resulting  temperature  profiles  showed  moderate  fits  by  all  methods  except  the  Full  Width 
method. 

The  strength  in  the  Fourier  Average  method  is  in  its  consistency  despite  range  of  data, 
or  data  elimination,  and  the  complete  rigor  in  which  the  parameters  are  calculated. 
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Appendix  A 

Fourier  Transform  of  a  Rectangular  pulse  and  Comparison  with  TRS  data 


The  Fourier  transform  of  a  time  dependent  equation  into  the  frequency  domain  is  defined 
as  such 

/+c© 

(A  - 1) 

►oo 

Given  a  discontinuous  function  described  as 


F(t)  —0,t  <  tt 


F{t)  =  A,U  <t  <tf 

Rt) 

A  _ 

T+- 

II 

o 

c*. 

IA 

Tl 

Tf 

Time 


then  integration  yields; 

(A-  2) 

were  i  =  >/— T  The  transform  evaluated  at  /  =  0  is  such; 

Jim  FU)  =  2ttA(1  -  l)/0 

(A -3) 

which  is  undefined.  Let 

gU)  =  i2xA(e~12irF‘  - 

(A -4) 

h(j >  =  2  nf 

(A -5) 

Then  differentiation  with  respect  to  /  yields; 

g'U)  =  2TrA(-Ue-i2*ft' 

(A -6) 

h'^ji  =  2n 

(A  -  7) 

Using  L’Hospital’s  rule  for  undefined  limits, 

Jim  =  A^tf  ~  ^ 

(A -8) 

Therefore, 

F{  o)  =  Apw 

(A -9) 

and  can  also  be  thought  of  as  the  area  under  the  curve  of  the  defined  function.  Using  Euler 
transformation, 


F(l)  = 


2  irf 


sin(2nfti)  —  sin(2vftj)  +  ilcosfixfti)  —  cos(2irftf)}  (  -  10) 
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and  defining  /o  such  that 

FUo)  =  0  (A  -11) 

This  can  only  happen  when  the  real  and  the  imaginary  parts  are  zero.  The  magnitude  of 
F(j )  is  defined  as; 

\\F\\  =  - — - Jreal 2  4-  imaginary 2 
27 r  f  y 


and  at  /  =  /0, 


only  if 


Ill’ll  =  2^JV2  ~  2cos[2irf{ti  -  tj)\ 

ll^ll  =  o 


2vf0(tj-ti)  =  2rn  n  =  1,2,3, ... 
if  tj  —  ti  is  the  pulse  width,  pw,  then  by  A-13, 


(A -12) 


(A-13) 


pu,  =  7o 


(>1-15) 


Thus  for  the  Fourier  Averaging  technique,  one  must  evaluate  the  fourier  transform  with 
/  =  0,  which  yields  the  fluence  of  the  Thermal  Radiation  Simulator  data.  Then  one  must 
find  the  minimum  value  of  the  magnitude  of  the  transform  and  determine  what  the  value 
for  /0  is.  This  yields  the  pulse  length,  which  leads  to  the  average  flux. 


One  could  take  this  type  of  analysis  farther  and  attempt  to  evaluate  the  starting  point  of 
the  TRS  record,  but  that  meaning  is  useless  in  this  aspect.  One  need  only  find  the  frequency 
in  which  the  real  or  the  imaginary  part  of  the  transform  is  zero,  then  compare  the  result  to 
a  rectangular  pulse  starting  at  t  =  0.  The  solution  is  quickly  derived,  but  withheld  in  this 
presentation. 
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Appendix  B 

Algorithm  Method  for  finding  Pulse  Width,  Fourier  Averaging 


Finding  the  pulse  width  is  the  key  to  the  successful  characterization  of  the  thermal 
pulse.  In  the  method  of  Fourier  Averaging,  the  pulse  width  is  determined  by  finding  the 
first  minimum  of  the  magnitude  of  the  fourier  transform.  Since  the  TRS  data  is  not  a 
perfect  rectangular  pulse  as  supposed  in  appendix  A,  the  minimum  may  not  be  zero.  So  the 
algorithm  used  must  search  for  a  minimum  value  as  opposed  to  a  zero  value  increasing  the 
difficulty  of  the  task.  Many  root  searching  algorithms  determine  a  zero  point  if  the  values 
searched  are  within  a  minimum  tolerance.  The  assumption  is  that  the  values  eventually  pass 
through  zero,  and  the  value  found  in  the  tolerance  is  sufficient.  Unfortunately,  the  TRS  data 
transform  minimum  may  be  greater  than  the  tolerance,  and  the  algorithm  would  pass  by  in 
sear  h  of  a  better  match. 


The  problem  is  solved  by  knowing  the  area  where  the  magnitude  minimum  is  located 
also  represents  a  drastic  change  in  the  values  of  the  magnitudes  derivative  with  respect  to 
the  frequency.  The  derivative  of  A- 12, 


ail*1l 

df 


A 

27 r/ 


£  - ~r(*/  ~  U)sin[2rf{tf  -  f,)] 

V2  -  2c°sP-*f(‘l  -  <‘)ll 


(B-l) 


shows  a  discontinuity  at  /  =  1  /(tj  —  tt)  =  f0,  A-15  as  shown  in  Figure  14.  The  TRS  data 
tram  forms  may  not  be  discontinuous,  but  the  point  where  the  change  in  signs  occurs  can 
be  found.  Initially,  the  point  /<j  is  guessed  at  by  approximating  the  pulse  width  as  that  of 


Figure  14.  Magnitude  and  Derivative  of  Fourier  Transform  of  50  cal/ cm2  —  sec  by  1  sec 
pulse. 


the  Full  Width  method.  Then  a  second  point,  /d  is  picked  nearby,  within  a  small  6f.  The 
respective  ||.F||  are  found.  The  results  are  used  to  determine  the  approximate  derivative  at 
fo  by; 


df  ~  /o2  -  /o1 


(B-2) 


A  third  point,  /q  is  initially  found  by  using  the  reciprocal  of  the  pulse  width  determined 
by  ‘  he  Full  Width  at  Half  Maximum  d?.ta.  a  fourth  point,  is  found  near  by  a  small 
6f.  The  magnitudes  are  evaluated,  and  the  approximate  derivative  of  is  found.  Then  the 
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results  of  the  two  derivatives  are  used  to  interpolate  for  /q.  A  point  near  /q  is  picked, 
within  some  6f,  and  the  magnitudes  and  approximate  derivative  of  /q  are  determined.  If  the 
derivative  of  /q  is  not  within  the  specified  tolerance,  the  process  is  repeated  with  /q  being 
replaced  by  the  closest  of  /o  or  /q. 


ft  pas  FI 

(tae  Full  Width  Method  an  First  tru. ) 


Find  Megl  and  Hag2 

FI  -F2 

(Use  Full  Width  at  Half  fadu  on  first  try) 

R  =  F3  +  <f  k- 


Find  Hag3  and  Hag4 


QQ: 

F3  -  R 

_  m 

F5  =  F1- 


0H1  (  FI  -  F3  ) 

Oil  -M 


freTiT# 


Find  Hag5  end  Mag6 


FI  =  F5 


DfQ 


F5  -  re 


No 


Is  1061  <  Tolerance? — No  —  Is  CH5  Negative  7 


Yes 


Yes 


Figure  15.  Flow  Chart  for  determining  Pulse  Width 
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Appendix  C 

Solution  of  the  Heat  Conduction  Equation  given  an  Initial  Temperature 

Distribution. 


For  the  1  dimensional  heat  conduction  equation 


d2Q  1  de 

dx2  a  dt 

(C-l) 

with  the  following  boundary  conditions, 

x  =  0,  ^  =  0 
ox 

(C-2) 

r  dQ  A 

X  =  L'  a7  =  0 

(C-3) 

and  initial  condition, 

t  =  0,  0  =  }(x) 

(C-4) 

a  solution  for  C-l  is, 

0  =  e-o2i(A1sm(Ax)  +  A2Cos(Ax)) 

(C  -  5) 

-rj-  =  e~“2t(A1cos(Ax)  +  A2sin(Xx)) 

(C  —  6) 

substituting  boundary  condition  C-2  sets  A\  =  0. 

Boundary  condition  C-3  sets  A2sin{XL)  =  0.  A2  —  0  leads  to  the  trivial  solution.  For 
the  nontrivial  solution,  sin(XL)  =  0,  only  if 


A„  =  — ,  n  =  1,2,3,...  (C  -  7) 

Thus  the  solution  by  boundary  conditions  is 

0  =  (c-8) 


the  initial  condition,  C-4,  will  determine  the  solution  of  An.  In  short,  An  is  determined  by 
use  of  the  orthogonality  relationship  for  a  cosine  function  [6].  A  countably  infinite  set  of 
functions  is  considered  orthogonal  in  the  interval  a  <  x  <  b  if 


( x)gn(x)dx  =  0  when 


m  =  n 


(C  -  9) 


let  /w  be  some  arbitrary  function  within  this  interval.  /(r)  can  be  expressed  in  terms  of 


fix)  =  Cigi(x)  +  C2g2(x)  +  ...  +  Cngn(x)  +  ...  +  Cmgm{ x)  +  ...  =  f^Cn0n(x)  (C  -  10) 


j=i 


The  C's  are  constants  to  be  determined.  If  the  series  C-10  is  convergent,  and  integrable  after 
multiplication  with  one  of  the  functions,  <7*(x),  then  by  definition  of  C-9, 


/  f(x)gk{x)dx  =  Ck  [  g\{x)dx 

J  o  J  a 


(C- 11) 
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(C  —  12) 


Thus  the  solution  for  heat  conduction  in  the  slab  after  exposure  to  the  thermal  load,  where 
w  is  the  displaced  time,  w  —  t  —  pw. 
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Appendix  D 

BASIC  and  Fortran  IV  Programs 

PARAM.FOR:  This  Fortran  program  was  the  instrument  useu  to  evaluate  the  four 
characterization  methods.  Note  that  the  data  are  all  based  on  a  sampling  rate  of  200  Hertz, 
and  no  lecord  is  over  five  seconds  long. 


CHARACTER  TRACE*15 
CHARACTER  METH*20 
COMMON  /  FLUX  /  F(4000),A(3) 

PRINT  Input  the  TRS  record  to  be  Characterized:’ 

READ  TRACE 

OPEN(UNIT=30,ERR=300,FILE=TRACE,STATUS=’OLD’) 

L1=0 

FMAX=0 

NSTRT=0 

NEND=0 

N=0 

5  N=N+1 

READ(30,100,ERR=11,END=12)  F(N) 

100  F0RMAT(F12.5) 

IF(F(N).GT.3.0)  THEN 
L  1=1 

GO  TO  10 
ELSE 

IF(Ll.EQ.l)  GO  TO  10 

NSTRT=N 

END  IF 

10  IF(F(N).LE.FMAX)  GO  TO  20 
FMAX=F(N) 

NMAX=N 

20  IF(F(N).GT.3.0)  GO  TO  5 

IF(N.LT.300)  GO  TO  5 
IF(NEND.NE.O)  GO  TO  5 
NEND=N 
GO  TO  5 
12  NSTP=N 
L2=0 

DO  200  N=NSTRT,NEND 
IF(F(N).GT.FM \X*0.5)  THEN 
L2=l 

NHME=N 

ELSE 

IF(L2.EQ.l)  GO  TO  200 

NHMS=N 

END  IF 
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200  CONTINUE 
SR=200 

C  FULL  WIDTH  METHOD 
NDN=NEND-NSTRT 
DN=REAL(NDN) 

SUM=0 

DO  220  N=NSTRT,NEND 
SUM=SUM+F(N) 

220  CONTINUE 
PWFW=DN/SR 
FLUE1=SUM/SR 
AFFW=FLUE1/PWFW 
C  FULL  WIDTH  AT  HALF  MAXIMUM 
NDN=NHME-NHMS 
DN=REAL(NDN) 

SUM=0 

DO  320  N=NHMS,NHME 
SUM=SUM+F(N) 

320  CONTINUE 

FLUE2=SUM/SR 
PWFWHM=DN/SR 
AFFWHM=FLUE2/PWFWHM 
C  MOMENTS  MATCHING  METHOD 
DO  410  1=0,2 
SUM=0 

DO  420  N=NSTRT,NEND 
T=REAL((N-NSTRT)/SR) 
SUM=SUM+F(N)*T**I 
420  CONTINUE 
A(I)=SUM/SR 
410  CONTINUE 
TC=A(1)/A(0) 

PWMMM=2*SQRT(3*(A(2)/A(0)-TC**2)) 

FLUE3=A(0) 

AFMMM=FLUE3/PWMMM 
C  FOURIER  AVERAGING 
PI=3. 1415927 
DF=0.001 
F5=1/PWFW 
450  Fl=F5 

F2=F1+DF 

CALL  SIMP(F1,NSTRT,NEND,AG1) 
CALL  SIMP(F2,NSTRT,NEND,AG2) 
DM1=(AG1-AG2)/(F1-F2) 
F3=1/PWFWHM 
455  F4=F3+DF 
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CALL  SIMP(F3,NSTRT,NEND,AG3) 

CALL  SIMP(F4,NSTRT,NEND,AG4) 

DM3=(AG3-AG4)/(F3-F4) 

460  F5=F1-DM1*(F1-F3)/(DM1-DM3) 

F6=F5+DF 

CALL  SIMP(F5,NSTRT,NEND,AG5) 

CALL  SIMP(F6,NSTRT,NEND,AG6) 

DM5=(  AG5-AG6)/ (F5-F6) 

IF(DM5.LT.-0.1)  THEN 

F1=F5 

F2=F1+DF 

CALL  SIMP(F1,NSTRT,NEND,AG1) 

CALL  SIMP(F2,NSTRT,NEND,AG2) 

DM1=(AG1-AG2)/(F1-F2) 

GO  TO  460 
ELSE 

IF(ABS(DM5).LT.0.05)  GO  TO  470 
F3=F5 
GO  TO  455 
END  IF 

470  PWFA=1/F5 

CALL  SIMP(0,NSTRT,NEND,FLUE4) 

AFFA=FLUE4/PWFA 

PRINT  *, ’METHOD  FLUENCE  PULSE  WIDTH  AVERAGE  FLUX  ’ 

METH=’FULL  WIDTH’ 

WRITE(*,500)  METH,FLUE1,PWFW,AFFW 
METH=’FULL  WIDTH  HALF  MAX’ 

WRITE(*,500)  METH,FLUE2,PWFWHM,AFFWHM 
METH=’MOMENTS  MATCHING’ 

WRITE(*,500)  METH,FLUE3,PWMMM,AFMMM 
METH=’FOURIER  AVARAGING’ 

WRITE(*,500)  METH,FLUE4,PWFA,AFFA 
500  FORMAT(1X,A20,1X,F12.6,1X,F12.6,1X,F12.6) 

STOP 

300  PRINT  *, ’ERROR  IN  THE  OPEN  STATEMENT  INPUT  FILE’ 

STOP 

11  PRINT  ’ERROR  READING  FILE’ 

STOP 

END 

SUBROUTINE  SIMP(FR,NSTRT,NEND,AG) 

COMMON  /  FLUX  /  F(4000),M(3) 

PI=3. 1415927 
SR=200 

THETA=FR*2*PI/SR 

SUMR=0 
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SUMI=0 

DO  600  N=NSTRT,NEND 
D=REAL(NSTRT-N+1) 
SUMR=SUMR+F(N)*COS(THETA*D) 
SUMI=SUMI+F(N)*SIN(THETA*D) 
600  CONTINUE 
RE=SUMR/SR 
RI=SUMI/SR 
AG=SQRT(RE**2+RI**2) 

RETURN 

END 


An  example  of  the  output  is  shown  below.  The  Data  input  was  Trace  3. 
Input  the  TRS  record  to  be  Characterized: 


TRACE3 

METHOD  FLUENCE 

FULL  WIDTH  10.659768 

FULL  WIDTH  HALF  MAX  9.182742 

MOMENTS  MATCHING  10.659768 

FOURIER  AVARAGING  10.659768 


PULSE  WIDTH  AVERAGE  FLUX 
0.700000  15.22824 

0.465000  19.747833 

0.576070  18.504309 

0.536578  19.866194 


MULTIPLT.BAS:  This  BASIC  program  was  used  to  generate  the  temperature  profiles 
using  series  of  artificially  generated  thermal  record.  The  output  was  read  to  an  ASCII  file 
called  "OUT?”  and  plotted. 


10  CLS 

20  DIM  T(100) 

30  KEY  OFF 
40  FOR  J  =  13  TO  24 
50  F$=”A:PULSE” 

60  G$=”D:OUT” 

70  IF  INT(J/10)=1  THEN  F$=”A:PULSEr 
80  IF  INT(J/10)=1  THEN  G$=”D:OUTr 
90  IF  INT(J/10)=2  THEN  F$=”A:PULSE2” 

100  IF  INT(J/10)=2  THEN  G$=,1D:OUT2,, 

110  J$=STR$(J) 

120  J$=RIGHT$(J$,1) 

130  FILE$=F$+J$ 

140  FILE2$=G$-f  J$ 

150  DUR=5 
160  SR=200 
170  INC  =  1 
180  N  =  SR  *  DUR 
190  R  =  1 

200  REM  **********  initialize  the  temperature  profile  ************* 
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210  FOR  K  =  1  TO  100 
220  T(K)  =  300 
230  NEXT  K 

240  REM  ***********  set  constants  ******************************** 

250  DT  =  1/200 

260  TI  =  0 

270  DX  =  .0015 

280  KO  =  204 

290  CP  =  896 

300  RHO  =  2707 

310  ALPHA  =  KO/(RHO*CP) 

320  A  =  .0001 

330  THETA  =  DT*ALPHA/DX2 
340  PHI  =  DX*41S60/I\O 

350  REM  **************  start  computation  **************************** 
360  OPEN  ”r,l, FILES 
370  OPEN  ”0”,2,FILE2$ 

380  PRINT” RUNNING  OUT  OF  FILE  ”;FILE$ 

390  PRINT#2,” TEMP  1  ” TEMP2” TEMP3” TEMP4” TEMP5” TIME” 

400  FOR  M  =  1  TO  N 

410  TI  =  TI  +  DT 

420  IF  EOF(l)  THEN  GOiD  460 

430  INPUT#1,Q 

440  IF  Q<2  GOTO  420 

450  GOTO  470 

460  Q  =  0 

470  T(l)  =  T(l)  +  THETA*(T(2)-T(1)+Q*PHI) 

480  FOR  K  =  2  TO  99 

490  T(K)  =  T(K)  +  THETA* (T(K-fl)  -  2*T(K)  +  T(K-l)) 

500  NEXT  K 

510  T(100)  =  T(100)  +  THETA*(T(99)-T(100)) 

520  IF  TI  <  R*.01  THEN  GOTO  550 
530  R  =  R  +  1 

540  PRINT#2,T(1),T(25),T(50),T(75),T(100),TI 
550  NEXT  M 
560  CLOSE#  1 
570  CLOSE#2 

580  PRINT  "End  Run.  Output  in  file  ”;FILE2$ 

590  CLS 
600  NEXT  J 
610  STOP 
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No  of 
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1 

Office  of  the  Secretary  of  Defense 
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Director,  Live  Fire  Testing 

ATTN:  James  F.  O’Bryon 

Washington,  DC  20301-3110 

1 

Director 

US  Army  Aviation  Research 
and  Technology  Activity 

ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 

2 

Administrator 

Defense  Technical  Info  Center 

Moffett  Field,  CA  94035-1000 

1 

ATTN:  DTIC-DDA 

Cameron  Station 

Alexandria,  VA  22304-6145 

HQDA  (SARD-TR) 

1 

Commander 

US  Army  Missile  Command 

ATTN:  AMSMI-RD-CS-R  (DOC) 
Redstone  Arsenal,  AL  35898-5010 

WASH  DC  20310-0001 

1 

Commander 

US  Army  Tank-Automotive  Command 

1 

Commander 

US  Army  Materiel  Command 

ATTN:  AMCDRA-ST 

ATTN:  AMSTA-TSL  (Technical  Library) 
Warren,  MI  48397-5000 

5001  Eisenhower  Avenue 

Alexandria,  VA  22333-0001 

1 

Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATAA-SL 

1 

Commander 

US  Army  Laboratory  Command 

White  Sands  Missile  Range,  NM  88002-5! 
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2 

Commander 
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1 

ATTN:  SMCAR-TDC 
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1 
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2 
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1 
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1 
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3 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3148 _  Date  Qf  Rcport  SEPTEMBER  1990 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  rcport  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


5.  Has  the  information  in  this  rcport  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Now  or  Corrcci 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


